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We study production of two pairs of jets in hadron-hadron collisions in view of extracting contri- 
bution of double hard interactions of three and four partons (3 — > 4, 4 — > 4). Such interactions, in 
spite of being power suppressed at the level of the total cross section, become comparable with the 
standard hard collisions of two partons, 2 — » 4, in the back-to-back kinematics when the transverse 
momentum imbalances of two pairing jets are relatively small. 

We express differential and total cross sections for two-dijet production in double parton collisions 
through the generalized two-parton distributions, 2GPDS [T], that contain large-distance two-parton 
correlations of non-perturbative origin as well as small-distance correlations due to parton evolution. 
We find that these large- and small-distance correlations participate in different manner in 4-jet 
production, and treat them in the leading logarithmic approximation of pQCD that resums collinear 
logarithms in all orders. 

A special emphasis is given to 3 — > 4 double hard interaction processes that occur as an interplay 
between large- and short-distance parton correlations and were not taken into consideration by 
approaches inspired by the parton model picture. We demonstrate that the 3 — » 4 mechanism, being 
of the same order in q s as the 4 — > 4 process, turns out to be geometrically enhanced compared to 
the latter and should contribute significantly to 4-jet production. 

The framework developed here takes into systematic consideration perturbative Q 2 evolution of 
2GPDS. It can be used as a basis for future analysis of NLO corrections to multi-parton interactions 
(MPI) at LHC and Tevatron colliders, in particular for improving evaluation of QCD backgrounds 
to new physics searches. 
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I. INTRODUCTION 

Understanding the rates and the structure of multi- 
jet production in hadron-hadron collisions is of primary 
importance for new physics searches. 

Production of high transverse momentum jets is a hard 
process which implies a head-on collision of QCD par- 
tons — quarks and/or gluons — from the small-distance 
wave functions of initial hadrons. Cross section of a 
hard collision is small compared with the size of hadron, 
a oc Q~ 2 ^hadr) wu; h Q 2 the scale related to trans- 
verse momenta of the produced jets, Q 2 ~ j\. Therefore, 
typically it is two partons that experience a hard colli- 
sion in a given event. A large angle scattering of these 
two partons produces two (or more) final state partons 
that manifest themselves as hadron large transverse mo- 
mentum jets. At the same time, one cannot exclude a 
possibility that more than one pair of partons happen to 
collide in a given event, giving rise to a multi-jet event. 
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A possibility of a double hard collision becomes more im- 
portant with increase of the energy of the collision where 
scattering off small x partons which have much higher 
densities becomes possible. 

In recent years multiparton collisions have attracted 
close attention. Following the pioneering work of Refs. [2, 
3 , a large number of related theoretical papers appeared 
j4j-[8] , based on the parton model and geometrical picture 
in the impact parameter space. More recently, this topic 
has been intensively discussed in view of the LHC pro- 
gram |TU] . Monte Carlo event generators that produce 
multiple parton collisions are being developed pTTTH3| : 
theoretical papers exploring properties of double parton 
distributions and discussing their QCD evolution have 
appeared dHHUl]. 

In our view, however, important elements of QCD 
that are necessary for theoretical understanding of the 
multiple hard interactions issue have not yet been prop- 
erly taken into account by above-mentioned intuitive ap- 
proaches. 

The problem is, sort of, educational: both the proba- 
bilistic picture, the MC generator technology is based 
upon, and the familiar Feynman diagram technique, 
when used in the momentum space, prove to be in- 
adequate for careful analysis and understanding of the 
physics of multiple collisions. 

From experience gained by treating standard (single) 
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hard processes, one became used to a motto that a large 
momentum transfer scale Q 2 ensures the dominance of 
small distances, r 2 ~ Q~ 2 , in a process under considera- 
tion. With the multiple collisions under focus, however, 
one has to distinguish two space-time scales: that of local- 
ization of the parton participating in a hard interaction, 
Ar 2 ~ Q~ 2 , and that of transverse separation, Ap, be- 
tween the two hard collision vertices. The latter can be 
large, of the order of the hadron size, even for large Q 2 . 

In order to be able to trace the relative distance be- 
tween the partons, one has to use the mixed longitudinal 
momentum-impact parameter representation which, in 
the momentum language, reduces to introduction of a 
mismatch between the transverse momentum of the par- 
ton in the amplitude and that of the same parton in the 
amplitude conjugated. 

Another unusual feature of the multiple collision anal- 
ysis that may look confusing at the first sight is the fact 
that — even at the tree level — the amplitude describing 
the double hard interaction process contains additional 
integrations over longitudinal momentum components; 
more precisely — over the difference of the (large) light- 
cone momentum components of the two partons originat- 
ing from the same incident hadron (see Section IV A[ ). 



In the previous short publication 1] we have considered 
production of two pairs of nearly back-to-back jets result- 
ing from simultaneous hard collisions of two partons from 
the wave function of one incident hadron with two par- 
tons from the other hadron ( "four-to-four" processes). 
As we have shown, this necessitates introduction of a 
new object — a generalized double parton distribution, 
2GPD, that depends on a new transverse momentum pa- 
rameter A conjugate to the relative distance between the 
two partons in the hadron wave function. Generalized 
double parton distributions provide a natural framework 
for incorporating longitudinal and transverse correlations 
between partons in the hadron wave function at the Qo 
scale, and for tracing the perturbative Q 2 -evolution of 
the correlations. 

The corresponding 4-jet cross section can be expressed 
in terms of 2GPD's as follows 



da(xi, X2, x 3 , £4) da 13 da 2A 
dti dt 2 dti dt 2 

d 2 A 



(1) 



(2tt): 



D a (xi,x 2 ; A) D b {x 3 , x 4 ; -A). 



The factor on the second line has dimension of inverse 
area: 



d 2 A D(x 1 ,x 2 ,A)D(x 3 ,x i ,-A 
nR 2 nt ~ J (2n) 2 D(x 1 )D(x 2 )D(x 3 )D(x i ) 



, (2) 



where D(xi) are the corresponding one-parton distribu- 
tions. The ratio of the product of two single-inclusive 
cross sections and the double-inclusive cross section 



(irR 2 nt ) is often referred to in the literature as "an ef- 
fective cross section" <j c g. We prefer, however, not to 
look at this quantity as a cross section, since it reflects 
transversal area of parton overlap as well as longitudinal 
correlations of the partons. At the same time, it has little 
to do with the measure of the strength of the interaction, 
which is what "cross section" represents. 

In a two-parton collision, scattered partons form two 
nearly back-to-back jets, while additional jets (should 
there be any) tend to be softer and to align with the di- 
rections of initial and final partons, because of collincar 
enhancements due to radiative nature of secondary par- 
tons. Such will be typical characteristics of a J^-jet event, 
in particular. On the other hand, four jets produced as 
a result of a double hard collision of two parton pairs 
would, on the contrary, form two pairs of nearly back-to- 
back jets. This kinematical preference is in stark contrast 
with "hedgehog-like" configurations of four jets stem- 
ming from a single collision and can be used in order 
to single out double hard collisions experimentally. 

Such experimental studies were recently carried out by 
the CDF and DO collaborations who have studied produc- 
tion of three jets + photon [Sl fTMTS] . The analysis of the 
data performed in [7] using information about general- 
ized parton distribution (GPDs) obtained from the study 
of hard exclusive processes at HERA has found that ob- 
served 3 jet + 7 rates were a factor > 2 higher than the 
expected rates based on a naive model that neglected 
correlations between partons in the transverse plane. 

The use of 2GPD allows one to incorporate such cor- 
relations and predict their Q 2 evolution. 

On the theory side, the back-to-back enhancement has 
been discussed, at tree level, in a number of studies 
of various channels (see, for example, discussion of the 
2 jets+66 in [19] and references therein). 

In the present paper we study perturbative radiative 
effects in the differential 4-jet distribution in the back-to- 
back kinematics and derive the expression for the corre- 
sponding cross section in the leading logarithmic collinear 
approximation. It takes into account QCD evolution of 
the generalized double parton distributions as well as ef- 
fects due to multiple soft gluon radiation, and turns out 
to be a direct generalization of the known "DDT formula" 
for back-to-back production of two large transverse mo- 
mentum particles in hadron collisions [20] . 

We also discuss and treat new specific correlations be- 
tween transverse momenta of jets due to 3-parton in- 
teractions producing 4 jets, "three-to-four" . Such pro- 
cesses are induced by perturbative splitting of a parton 
from one of the hadrons, the offspring of which enter 
double hard collision with two partons from the wave 
function of the second hadron. The hard scale of this 
parton splitting is determined by transverse momentum 
imbalances of pairs of jets, S\ 3 , 5 2 a, and exhibits specific 
collincar enhancement in the kinematical region where 
two jet imbalances practically compensate one another, 
5 12 = (S 13 + 8 2i ) 2 « S 2 13 ~ S 2 24 . 
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Consistently taking into account three-to-four parton 
process solves a longstanding problem of double counting 
in treating multi-parton interactions. 

Discussion of the 2-parton distribution has a long his- 
tory. It is commonly defined in the momentum space as 
a 2-particle inclusive quantity depending on two parton 
momenta, see [5TJ [35] . Being related to (the imaginary 
part of) a certain forward scattering amplitude, it there- 
fore disregards impact parameter space geometry of the 
interaction. Exploring properties of 2-parton distribu- 
tions so defined, an approach to the study of the multiple 
jet production has been recently suggested in Ref. [23] . 
The reason why this approach has faced difficulties, [2~4"] . 
and did not solve, in our view, the problem of system- 
atic pQCD analysis of 4-jet production is clear: it did 
not incorporate effects due to variations of the transverse 
separation between the partons — information encoded 
by 2GPD's but not by the 2-parton momentum distribu- 
tions. 

The 2GPD's were recently used in Ref. [25] for intuitive 
description of the total 4-jet production cross section. 
However, the differential distributions were not discussed 
in that paper, and not all relevant pQCD contributions 
were included, so that our results arc different from the 
ones obtained in [25] . 

The paper is organized as follows. 

In Section [TT] we recall the main ingredients of the 
perturbative analysis based on selection of maximally 
collincar enhanced contributions in all orders. In Section 
|III| we present the evolution equation for generalized two- 
parton distributions. Section |IV] is devoted to the per- 
turbative analysis of small-distance correlations between 
partons. The main result of the paper — the differential 
distribution of 4-jet production in the back-to-back kine- 
matics — is formulated in Section |Vj and the total cross 
section of two-dijet production is described in Sect ion|VI| 
Conclusions and outlook are presented in Section VII 



introduce transverse momentum imbalance: ji± +js± — 
6 is which constitutes another hard scale: Q\ 3> 5f 3 3> 
Aq CD . For production of four jets in the back-to-back 
kinematics, this gives four different hard scales. As we 
shall see below, in the 3-partons collisions producing four 
jets yet another scale enters the game: 5' 2 3> Aq CD with 

5' = 5i3 + 824 — the total transverse momentum of the 
4-jet ensemble. 

In what follows we consider transverse momenta of all 
four jets to be of the same order, Q 2 = jf x ~ j 2 ± ~ 
j| ± ~ j 2 j_. This is not necessary but helps to avoid 
complications in the hierarchy of relevant scales. 

Finally, let us mention that in what follows it will be 
tacitly implied that fixing these scales — from the largest 
one, Q 2 , down to smaller ones, S and §' ', — is not compro- 
mised by uncertainties in determination of the transverse 
momenta of the jets. 



B. Back-to-back kinematics 

The basic 2-jet production cross section scales, asymp- 
totically, as 



da { 



2^,2) 



dt 



(3) 



According to Q, production of four jets in simultaneous 
hard collisions of four partons yields 



da^ 
dtidt 2 



cx R- 



R 2 Q 8 



(4a) 



with R 2 = 1/ (A 2 ) the characteristic distance between 
the two partons in the hadron wave function. At large 
Q 2 this cross section is parametrically smaller than that 
for production of four well separated jets with transverse 
momenta jfj_ ~ Q 2 in a 2-parton collision: 



II. PERTURBATIVE ANALYSIS 
A. Hard scales 

The perturbative approach implies that all hard- 
ness (transverse momentum) scales that characterize the 
problem are comfortably larger than the intrinsic QCD 
scale Aqcd: Q\ 3> Aqcd- The process under considera- 
tion may have up to five hard scales involved. 

Indeed, in the leading order in a s , large transverse mo- 
mentum partons are produced in pairs and have nearly 
opposite transverse momenta, setting the hard scale 
Q\ = jf x ~ j 2 ± . Within the parton model framework 
(neglecting finite smearing due to intrinsic transverse mo- 
menta of incident partons), one has da cx <5(j'i_l + jax)- 
Secondary QCD processes — evolution of initial parton 
distributions and accompanying soft gluon radiation — 



riq( 2 ^ 4 > 
diidi 2 



Q 6 



(4b) 



(with transverse momenta of two out of four jets being 
integrated over). 



Qualitatively, the production mechanism (4a I can be 
labelled a "higher twist effect". Nevertheless, it may 
turn out to be essential — comparable with the "lead- 
ing twist" 2 — > 4, Eq. (4b) — if one looks at specific 



kinematics of the 4-jet ensemble. 

Let z be the direction of colliding hadron momenta. 
Imagine that we are triggering on two jets moving along 
the x and y axes in the transverse plane, and look for 
two accompanying jets inside some solid angles Afi <C 
47r around the —x and — y directions. The production 



mechanism (4b I docs not populate this region: the higher 



order 2 — > 4 QCD matrix element is enhanced when two 
final state partons become quasi- collinear but is perfectly 
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smooth in the back-to-back kinematics. Therefore, its 
contribution will be suppressed, 



(An 



1 



R 2 Q 2 ' 



D. Double Logarithmic parton form factors 

In the leading order in a s , it suffices to have just one 
parton present with k± = —813 in order to assure J13 7^ 0. 
At the same time, inclusive production of accompanying 
partons with transverse momenta k± turns out to be sup- 



contrary to the 4 ^ 4 production mechanism f4a| which presged in & brQad interyal « fc 2 <<; Q 2 j as long ag 



is concentrated in this very kinematical region 



C. Collinear approximation 

The differential 4-jet production cross section possesses 
two collinear enhancements. Depending on the kinemat- 
ics of the jets, they are, symbolically, 



dcr (4 ^ 4) oc W~72" d 2 j3±d 2 j 4± ■ dY,, 



(5a) 



§12 §2 



S' 2 « S 2 « Q 2 , 6 2 = 6j 3 ~ S 2 4 . (5b) 



Here dT, = d£(ii,i2) is the cross section integrated 
over the transverse momenta of the "backward" jets 3 
and 4. The integrated cross section dS contains the 
squared matrix element of the four-parton production 
and is of the order of a 4 , cf. Eq. (4a). At the Born 



level, the jets in pairs are exactly back-to-back, so that 
da^ 4 ^^ cx (5 2 (<5i3)<5 2 (<524) in (Sal. To have a non-zero 
value of the transverse momentum imbalance, one has 
to have additional large transverse momentum parton(s) 
produced. 

In the second important contribution to the cross sec- 
tion, Eq. (5b), one power of the coupling emerges from 



the splitting of a parton from one of the incident hadrons 
into two, and the second power is due to production of 
an additional final state parton with k± = —5'. 

In both cases the smallness due to additional powers 
of the coupling is compensated by two broad (logarith- 
mic) integrations over transverse momentum imbalances 
as indicated in §5y. 



one wants to preserve the collinear enhancement factor 
<5 2 3 in the jet correlation (5a). 



This dynamical "veto" has two consequences. 

First of all, it results in reduction of the hardness scale 
of the parton distributions from the natural scale Q 2 
(scale of the parton distributions in the integrated cross 
section) down to the observation-induced scale <5 2 3 <C Q 2 ■ 

Then, it introduces double logarithmic (DL) form fac- 
tors of participating initial state partons, since the trans- 
verse momentum of the jet pair can be compensated not 
only by a hard (energetic) parton from inside initial par- 
ton distributions but also by a soft gluon whose radiation 
did not affect inclusive parton distributions due to real- 
virtual cancellation. 

The presence of the DL form factors depending on the 
logarithm of a large ratio of scales, \n(Q 2 /S 2 A, is typical 
for the so-called "semi-inclusive" processes [201 ■ 

Production of massive lepton pairs in hadron collisions 
(the Drell-Yan process) is a classical example of a two- 
scale problem. Here enter form factors of colliding quarks 
that depend on the ratio of the invariant mass q 2 to the 
transverse momentum of the lepton pair, a s \n 2 {q 2 /q\), 
in the dominant kinematical region q\<SL q 2 : 



da 



dq 2 dq 2 . 



da to 
~dq^ 



Dl( Xl , q l)Dl( X2 ,ql)S 2 q (q 2 ,ql) .(6) 



dq 2 



S q is the double logarithmic QCD quark form factor. 

Sudakov quark and gluon form factors can be expressed 
via the exponent of the total probability of the parton 
decay in the range of virtualities (transverse momenta) 
between the two hard scales: 



S q (Q 2 ,n 2 ) = exp<| - j 2 f^^p fj k ' Q dzP«{z) \ , (7a) 



S g {Q 2 ,. 2 ) = cxp<|- f?^ 1 [ m dz[zPI{z)+n f PI{z)]\. (7b) 







Here P^iz) are the non-regularized one-loop DGLAP splitting functions (without the "+" prescription): 
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l + z A 



P°{z)=P*(l-z), 



Pl{z) = T R [z 2 + (l-z) 2 ] , P|(*) = C A 



l+z 4 +(l-z) 4 
z(l-z) 



(8) 



r 



the upper limit of z-integrals properly regularizes the soft 
gluon singularity, z — > 1 (in physical terms, it can be 
looked upon as a condition that the energy of a gluon 
should be larger than its transverse momentum, [20] ) . 

The case of hadron interactions producing large trans- 
verse momentum partons (instead of colorless objects like 
a Drell-Yan pair or an intermediate boson) is more in- 
volved since here the transverse momentum imbalance 
may be compensated by QCD radiation from the final 
state partons too. 

The azimuthal correlation between two nearly back- 
to-back large transverse momentum particles was consid- 
ered in [5D]. An analog of the "DDT formula" has been 
derived in the collinear approximation, which expression 
contained the product of four form factors, two initial 
parton distributions and two fragmentation functions. 



of the parton system, see [27] and references therein. 

In the present paper we concentrate on resummation 
of collinear enhanced DL and SL terms and avoid com- 
plications due to soft SL corrections. This means ig- 
noring color transfer effects in hard interactions. Thus, 
production of four jets with large transverse momenta 
j± ~ Q and pair imbalances will be equivalent, in our 
treatment, to production of two colorless Drell-Yan pairs 
with invariant masses O (Q) and transverse momenta ^13 
and 624- Generalization of the results of (27] to the case 
of double parton scattering seems straightforward and 
should be considered separately. 



III. GENERALIZED DOUBLE PARTON 
DISTRIBUTION 

A. Geometry of 2 GPD 



E. Single Logarithmic soft gluon effects 

The case when jets are being reconstructed in the final 
state is more complicated to analyze as it yields an an- 
swer depending on the jet finding algorithm. The prob- 
lem has been addressed by Banfi and Dasgupta in [27] 
where a smart way of defining the final state jets was 
formulated that permitted to write down a resummed 
QCD formula for soft gluon effects in "2 partons — » 2 
jets" cross sections. 

Collinear logarithms due to hard splittings of the fi- 
nal state partons do not pose a problem: such secondary 
partons populate the jets. Partons that appear as sep- 
arate out-of-jet radiation — and are relevant for trans- 
verse momentum imbalance compensation — have to be 
produced at sufficiently large angles with respect to the 
jet axis. This is the domain of large-angle gluon radi- 
ation. Production of soft gluons in-between jets is also 
logarithmically enhanced and induces single logarithmic 
(SL) corrections, [a s ln(Q 2 /<5 2 )] , that may also be sig- 
nificant and should be resummed in all orders. 

Contrary to collinear enhanced effects (that drive evo- 
lution of parton distributions and fragmentation func- 
tions and determine the Sudakov form factors), the large- 
angle gluon radiation cannot be attributed to one or an- 
other of the partons participating in the hard scattering. 
It is coherent and depends on the kinematics and color 
topology of the hard parton ensemble as a whole. As a 
result, resummation of these SL corrections becomes a 
matrix problem that involves tracing various color states 



The 2GPD in the expression for the multiparton pro- 
duction cross section has a meaning of a two body form 
factor when partons 1 and 2 receive transverse momenta 
A and —A leaving the hadron intact. Nonrelativistic 
analogue of this form factor is familiar from the double 
scattering amplitude in the momentum space representa- 
tion of the Glauber model, see e.g. [25]. Recall that [T] the 
scale A in 2GPD is conjugate to the relative transverse 
distance between the two partons in the 2GPD in the 
impact parameter representation considered in [2] [3] 114) . 

Two partons may originate from soft low-scale fluc- 
tuations inside the hadron; they can also emerge from 
a perturbative splitting of a common parent parton at 
relatively large momentum scales. It is clear that these 
two contributions to 2GPD will have essentially different 
dependence on the parameter A. 

The first contribution we will denote 

[2 p b a c (x 1 ,x 2 ;ql,ql,A), (9) 

with the subscript [ 2 ] stressing that here the partons b 
and c emerge from the no-perturbative wave function of 
the hadron a. It should decrease rapidly at scales larger 
than a natural scale of short-range parton correlation in 
a hadron (this scale may be slightly different for quarks 
and gluons and could in principle be significantly larger 
than the 1/rjv as there exists another non-perturbative 
scale of the chiral symmetry breaking which maybe as 
large as 700 MeV). 



6 



2GPD should rapidly decrease for A 2 > 1.5/ (r 2 g ) 



where (i 



\l/2 



is the transverse gluonic radius of the nu- 



cleon. In the mean field approximation when the corre- 
lations are neglected, it can be approximated by a fac- 
torized expression |T] 

[2 p b a ' c (x u x 2 ;qlqlA) = F g (A 2 ;x 1 ,q 2 1 )F g (A 2 ;x 2 ,q 2 2 ) 
x G b a {x u q 2 )G c a (x 2 ,q 2 ), (10a) 

where G are the single-parton distributions and the two- 
gluon form factor F can be parametrized as 



F g {A 2 



A 2 



m 2 (x) 



(10b) 



The parameter m 2 is of the order of 1 GeV 2 for x ~ 10~ 2 
and gradually drops with decrease of x and with increase 
of virtuality |29j . 



The second contribution we will denote 



[il 



pY{xx,x 2 ;ql,qi;A), 



(11) 



where the subscript ^ stands as a reminder of the fact 
that b and c originate from perturbative splitting of a sin- 
gle parton from the hadron wave function. This contri- 
bution is practically A-independent and should decrease 
with A much more slowly, due to logarithmic pQCD ef- 
fects. (A steep power falloff starts only when A 2 exceeds 
the relevant hard scale, q 2 .) 



B. On the geometrical enhancement of the 
interference effects due to pQCD correlations 



By total cross section in the present context we mean 
the back-to-back 4-jet cross section integrated over pair 
jet imbalances in the dominant logarithmic region Sik <C 



We start by noting that the product of two small- 
distance parton fluctuations, \ip x ^jD, does not con- 
tribute to the process we are interested in. Indeed, in 
this case the integral over A in Eq. (fTl) formally diverges 
and yields a hard scale (instead of in the numerator. 
This means a significant contribution to the cross section 



but not the one we are looking for. Below in Section IV B 
we will explicitly verify that a double hard collision of two 
parton pairs each of which originates from perturbative 
splitting, lacks the back-to-back enhancement. In fact, 
the product x [ijD corresponds to a one-loop cor- 
rection to the "leading twist" perturbative production of 
four jets in a hard collision of two partons ( "two-to-four" ) 
whose distribution is smooth in the back-to-back region 
and as such gets subtracted as background. 



Keeping this in mind, the back-to-back 4-jet produc- 
tion cross section is proportional to the inverse "interac- 
tions area" S described by the expression 



= J |^(p)D a (A) [2] A,(A) + [2] D (A) {1 p b (A) + {1 p a (A) [2 p b (A)\ 



(12) 



where indices a, b mark two interacting nucleons. This 
expression is somewhat symbolic; a careful analysis of 
the "interaction area" will be carried out below in Sec. [V] 



(see Eq. (|32j)). 

The first term in Eq. ( 12 ) we will refer to as a "four- 
to-four" process: two partons from the wave functions of 
the hadron a interact with two partons from the hadron 
b producing four jets. The second and the third terms in 
Eq. ( 12 ) describe hard collisions of one parton from one 
hadron with two partons from the second hadron. Until 
recently, these "three-to-four" processes were commonly 
ignored in the literature (see, however, [3S] ). At the 
same time, they turn out to be somewhat enhanced. 

Indeed, the contribution due to four-to-four processes 



to the geometrical factor Eq. ( 12 ) is given by 



£A_ 



F 2 (A 2 ) x F 2 (A 2 ) 



rn 



9 

7ir ' 



(13a) 



leads to fast convergence of the integral whose median is 
positioned at a value as low as A 2 w 0.1m 2 . 

The case of the three-to-four process is different. This 
process corresponds, as we explained above, to interac- 
tion of the offspring of the perturbative splitting of a 
parton from the wave function of one hadron, with two 
partons from the non-perturbative wave function of the 
second colliding hadron. On the side of m-D, the param- 
eter A enters "perturbative loop" due to parton splitting 
and as a result the dependence of m Don A turns out to 
be only logarithmic, that is, parametrically much slower 
than that of the non-perturbative form factor F g (A 2 ). 
Thus, the three-to-four contribution to the double inter- 
action cross section reduces to 

fd 2 A , , fd 2 A ,, A , 

J w^* 1 '* 25 A) ^ (A ^ m D \±*J (2^ F * (A ' 

(13b) 



), 



Fast decrease of the product of two squared form factors (where we have neglected the logarithmic A-dependence 
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of m-D). This corresponds to the fact that in the impact 
parameter space, the distance between partons coming 
from a perturbative splitting is much smaller than the 
hadron size, so that the answer is proportional to the den- 
sity of non-perturbative two-parton correlation at small 
distances — "in the origin" : 



dA 9/a9 . 

F 2 (A 2 ) 

(2tt) 2 g{ ' 
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3tt ' 



(13c) 



Comparison with the estimate (13a) shows that the con- 
tribution to the cross section of the "interference term" 
1 + 2 — » 4 is enhanced, relative to the 2 + 2 — > 4 process, 
by the factor 



7 



x 2 - 5. 



(14) 



(For the case of the Gaussian form factors this enhance- 
ment is 15% smaller — a factor of 4.) This estimate was 
obtained for the case when all four partons participating 
in the hard collisions are gluons. A detailed numerical 
study of the a^-dependence of the effective interaction 
area S will be presented in the paper under preparation. 

So we conclude that the three-to-four processes may 
provide a sizable contribution to the cross section even if 
they constitute a small correction to 2GPD. 



C. Perturbative QCD effects in 2 GPD 



Thus, we represent the generalized double parton dis- 
tribution 2GPD as a sum of two terms: 



D b a > c {x 1 ,x 2 ;qlq 2 ;&) = [2p b a ' c { Xl , x 2 ; q\, q\\ A) 

+ ll p b a ' c (x u x 2 ;qlq 2 2 ;A). (15) 



The term [ 2 ] D describes the distribution of two partons 
from the non-perturbative wave function of the hadron 
a that are independently evolved to large perturbative 
scales q\ and q\ according to the standard one-parton 
evolution equation. The perturbative evolution involves 
momentum scales much larger than the hadron wave 
function correlation scale (A 2 ) ~ Q\. Therefore, the evo- 
lution practically does not affect the A-dependence of the 
two-parton spectrum. This piece of the 2GPD acquires 
but a mild additional logarithmic dependence at the tail 
of the A-distribution in addition to a non-perturbative 



power falloff ( 10 1 



The integral QCD evolution equation for [ 2 ]D bc reads 



[2] D b a c ( x 1: x 2 ;qlql A) = S b (ql Q 2 min )S c (ql Q 2 min ) [2] D b f(x 1 , x 2 ; Q 2 , Q 2 ; A) 



Here Pf(z) are the non-regularized one- loop DGLAP 
splitting functions pi and Si — the double logarithmic 
Sudakov parton form factors defined in Eq. ([7]). 

The lower limit of the perturbative evolution in 
Eq. p|, 



max(Q 2 ,A 2 ) ~ Q 2 + A 2 



(17) 



is the only source of additional (logarithmic) A- 
dependence. It emerges when A 2 exceeds — and sub- 



stitutes — the starting non-perturbative scale Qq of the 
perturbative evolution. 



The second term in Eq. ( 15 ) 



represents the 



small-distance correlation between the two partons that 
emerge from a perturbative splitting of a common par- 
ent parton taken from the hadron wave function. It can 
be expressed in terms of standard inclusive single-parton 
distributions as follows 



ll] D b -(x 1 ,x 2 ;qlq 2 ;A)= £ \% lr Ag , 

a' .0' ,G ^mm 



2 /-)2\ 

<!.*) 

f dz ^h'\r'}/ \ ( 3>1 9 , 9 \ ( %2 9 , 9 \ 



P b a y\z) G\, (^;q 2 ,k 2 ) G c c , ( jr^-i&k? 
z{l-z) \zy J \{l-z)y 



The A-dependence of mD is very mild as it emerges 
solely from the lower limit of the logarithmic transverse 
momentum integration Q 2 



) 2 . . 

■ruin 



IV. ANALYSIS OF PERTURB ATIVE 
TWO-PARTON CORRELATIONS 



Let us analyze the lowest order interaction amplitude 
shown in Fig.[T]that produces a double hard collision and 
involves parton splitting. 




enter the denominator of the amplitude in terms of the 
Sudakov variables read 

k\ = x\(3s — k\ — —x 2 (3s — k\, 

with s — 2p a pb and k\ = (fcjj 2 > the square of the 
two-dimensional transverse momentum vector. 

A singular contribution we are looking for originates 
from the region (3 -C 1, so that precise form of the longi- 
tudinal smearing does not matter and the integral yields 



N 



dd 



2-KiN 1 



[x\(3s — k\ +ie)(— x 2 f3s — k\ +ie) (X1+X2) k\ ' 

The numerator of the amplitude is proportional to the 
first power of the transverse momentum k±. As a result, 
the squared amplitude (and thus the differential cross 
section) acquires the necessary factor 1/d 2 that enhances 
the back-to-back jet production. 



B. 



Now we should verify that the diagram of Fig. [2] where 
both incident partons split, and their offspring engage 
into double hard scattering, does not favor back-to-back 
jet kinematics. In other words, it does not lead to a small 
imbalance factor 1/S 2 in the differential cross section. 



FIG. 1: 3^4 

We express parton momenta fc, in terms of the Sudakov 
decomposition using the light-like vectors p a , pb along the 
incident hadron momenta: 




k\ = xip a + fipb + k±, k 3 ~ (i 3 - (3)p b ; 
k 2 = x 2 p a ~ PPb ~ k±, fe 4 ~ (i 4 + j3)p b ; 
k ± = S 12 = -<5 34 (<5' = 0); fc ~ (x\ + x 2 )p a . 

Here fc , k 3 and ki are momenta of incoming (real) par- 
tons, and fci and k 2 — virtual ones. Light-cone fractions 
Xi, i=l,..4, are determined by jet kinematics (invariant 
masses and rapidities of jet pair). The fraction j3 that 
measures the difference in longitudinal momenta of the 
two partons coming from the hadron b, is arbitrary. Fixed 
values of the parton momenta x% — (3 and X4 + (3 cor- 
respond to the plane wave description of the scattering 
process in which the longitudinal distance between the 
two scatterings is arbitrary. This description does not 
correspond to the physical picture of the process we are 
interested in. In order to ensure than the partons 3 and 
4 originate form the same hadron of finite size, we have 
to introduce an integration over (3 in the amplitude, in 
the region f3 = O (1). 

The Feynman amplitude contains the product of two 
virtual propagators. The virtualities k\ and the k\ that 



3 \ / 4 
5 " 

FIG. 2: 2^4 
Sudakov decomposition: 

ki = (xx-a)p a + f3p b + k' ± , k 3 = (x 3 -f3)p b + ap a - k±; 
k 2 = (x 2 +a)p a - [3p b - k' ± , ki = {x A +(3)p b - ap a + k ± ; 

k' ± -k ± =S 12 = -S 34 (S' = 0); 
k ~ (xi + x 2 )p a , k 5 ~ (x 3 + x 4 )p b . 

This is a loop diagram and it contains explicit integration 
over the loop momentum: 

f dad(3 f d 2 k ± d 2 k' ±!;2/r - r , . 
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To get an enhanced contribution we have to have parton 
virtualities that enter the denominator of the Feynman 
amplitude to be relatively small, of the order of S 2 <C Q 2 . 
This implies \a\, <C 1 in the essential integration re- 
gion. Adopting this approximation, we can simplify par- 
ton propagators and reduce the longitudinal momentum 
integrations to the product of two independent integrals: 



dp 



2ni (j3x\S — k\ + ie)(l3x2S + k 2 _ — it) 
da s 
2ni (ax 3 s~k 12 + ie)(ax4S + k' 2 — ie) 
i 1 



(xi + x 2 )(x 3 + x A )s k\ kf ' 

The remaining transverse momentum integration takes 
the form 



d 2 k ± V 



(20) 



Due to gauge invariance the numerator of the diagram — 
the "vertex factor" V — is linear in transverse momenta 
of the loop partons: V oc k±k±. Therefore, the integral 
( 20 1 produces no more than a logarithmic enhancement 



factor, \n(Q 2 / 5 2 ) , instead of the power back-to-back sin- 
gularity Q 2 /5 2 we were looking for. 

So, the diagram Fig. [2] with double parton splitting 
constitutes but a negligible loop correction to the usual 
"hedgehog" 4-jet kinematics typical for 2 — !• 4 QCD pro- 
cesses. The fact that this loop diagram does not produce 
a pole singularity in 5 2 could have been extracted, e.g., 
from numerical studies of double Z-boson production in 
two-parton collisions |30j and, more generally, of multi- 
leg parton amplitudes [31]. 

The logarithmic character of this correction has been 
recently confirmed by the systematic study of "box inte- 
grals" in gl]. 

The presence of the double parton splitting contri- 
bution of Fig. [2] is being treated in the literature as a 
source of potential problem of double counting (see, e.g., 
[IH1 HU H2] ) . The present paper solves this problem. 



C. 3 — > 4 with additional parton emission 

We have to return now to the 3 — ¥ 4 process and exam- 
ine the possibility of producing an additional parton, in 
collinear enhanced manner, in order to lift off the Born 
level kinematical constraint 5' = 0. 





(a) 



(b) 



FIG. 3: Three-to-four amplitudes with extra emission from inside the splitting fork 



Consider the diagram of Fig. [3^,. The momenta of We have three virtual propagators subject to integration 
quasi-real colliding partons are over /3: 



S' 2 

k ~ (x 1 +x 2 +a)p a ; k 3 ~ (x 3 -/3)p b , k 4 ~ (x 4 +/3H )p b , 

as 



ki = xip a + /3p b + 5n, (21a) 



and the radiated on-mass-shell parton carries momentum , / _ S' 2 , 

k 2 = X2Pa - [P+—)pb + $24, (21b) 



<5'2 V as , 

aPa + ^s~ Ph ~ 5 '' 5 ' = ^ + 5 ^- i = (x 2 + a)p a -j3p b -6 13 . (21c) 
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Closing the contour around the pole k\ 
obtain (3 s = S\ 3 jx\ and 



- e = (x 2 + a)/3s + 5 

-hi = x 2 — )s 
\ as J 



S 2 



,,^±^, (22a) 

2 + ^^ 22b) 



Taken tog ether with the residue of the /3-integration, 
Eq. (22b) produces the universal factor present in 
all the amplitudes considered, (including the diagrams 
with parton emission off the external lines, see Fig. [4] 
below): 



P- 



xi<5 24 



2 X\X2 r.,2 
X20 13 H . 



(23) 



the parton 1 collides with 3, while the parton £ keeps 
evolving and scatters off 4 with a much larger momen- 
tum transfer 5 2 4- Evolution of the parton £ in between 
these two scales is the origin of probable (logarithmically 
enhanced) production of additional parton(s). 

Analogously, the diagram Fig. [3]d with a parton pro- 
duced off the virtual line 1 contributes in the complemen- 
tary kinematical region 



4 



<. 6 



13 



(24b) 



Full perturbative analysis of the production of a parton 
from inside the "splitting fork", Fig. |3j together with 
emission off the incoming line "0" (to be treated below 
in Sec. VB) is sketched in the Appendix. 



We observe that both propagators (22) are enhanced in 
the back-to-back kinematics. The amplitude of Fig. [3^ 
gives a double collinear enhanced contribution to the 
cross section in the region 



'13 



(24a) 



The inequality (24a) corresponds to the following phys- 
ical picture. An incident parton fco splits into k\ and £ 
early, at time O (y/s/df 3 ) corresponding to some com- 
paratively low perturbative scale <5i3 . At this time scale 



V. DIFFERENTIAL DISTRIBUTION 

A. 2 + 2^4 

Now that we know the structure of the perturbative 
corrections to 2 GPD, Eqs. ( 15 )-( 18 ), we are in a position 
to write down the generalization of the DDT formula ^ 
for (the first contribution to) the differential cross section 
of 4-jet production in nearly back-to-back kinematics. It 
reads 



da<- 4 ^ da Dart d d 



d 2 5 13 d 2 S 24 dhdh dS 2 13 dS 2 4 



Tj-2 uu _ '"V'i'i v v J ...r,l,2f„. ~,-.JC2 x 2 ^ „ . ,r>3,4/_ „ . . x 2 X 2 



[2pa {xi,x 2 ;S l3 ,5 24 ) X [2] ZV (x 3l X 4 ] S 13 , S 24 ) 



Si (Q 2 ,S 2 13 ) S 3 {Q\5 2 13 ) x S 2 (Q 2 , S 2 24 ) S 4 (Q 2 , <5 2 4 ) . (25) 



Here dcr par t is the cross section of double hard parton 
scattering, and Si stand for Sudakov form factors of four 
participating partons. Sum over parton species and con- 
volution over A as in Eq. ([I]) is implied. 

Taking derivative over the scale S 2 of the function de- 
pending on a s \ogS 2 , produces the factor a s /6 2 . Differ- 
entiating the Sudakov form factor of a given parton de- 
scribes the situation when the jet imbalance is compen- 
sated by radiation of a soft gluon off this parton. Dif- 
ferentiation of the parton distribution corresponds to the 
situation when a hard parton takes the recoil. 



B. 1 + 2^4 



The differential transverse momentum imbalance dis- 
tribution due to the cross-terms [ 2 iD x contains two 
pieces. 

1. Two compensating partons 



The first one has the same structure as Eq. ( 25 ) 



da\ 



(3-Hl) 



do 



part 



d d 



lipa 2 (x 1 ,x 2 ;S'f 3 ,S^ 4 ) ■ m D 3 b ' 4 {x 3 ,x 4 ;5'f 3 ,S^ 4 ) 



d 2 S 13 d 2 5 24 di l dt 2 dS 2 3 dS 2 4 

x Si (Q 2 ,5 2 13 ) S 3 (Q 2 , 5 2 3 ) ■ S 2 (Q 2 , <5 2 4 ) S 4 (Q 2 , 5 2 4 ) 



(26) 
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Sum over parton species and convolution over A is im- 
plied as above in Eq. ( |25] ). 

As above, taking the derivatives in Eq. (26) corre- 



sponds to fixing transverse momenta of two final state 
partons that compensate jet pair imbalances: — S13 and 

-^24- 

Consider now the correlation term m D of Eq. ( 26 1 . If 



we apply the derivatives to the parton distributions D in 



the integrand of Eq. (18 1 for the correlation term m-D, 
this contribution will correspond to production of two 
momentum compensating quanta in the course of evolu- 
tion of the system of two partons with account of small- 
distance perturbative correlation between them. In this 
case the scale of the core parton splitting stays smaller 
than the two external scales #13, #24) and is being inte- 
grated over. 



2. One compensating parton 



The correlation term ( 26 1 contains an additional op- 



tion. Namely, instead of creating intermediate state par- 
tons that keep evolving up to external scales, the pertur- 
bative splitting may produce that very parton that gets 
engaged in the hard scattering. This possibility is also 
contained in Eq. ( 26 1 : it corresponds to the differentia- 



tion of the upper limit of the virtuality integral in Eq. ( 18 ) 
over the smaller of the two imbalances, min{<5^ 3 , S^}- 
One of the two parton distribution functions in the in- 
tegrand then collapses to 5(1 — x/y). Taking the second 
derivative over the larger imbalance of the second par- 
ton distribution produces the contribution described by 



the diagram of Fig. [3] that we have discussed above in 
Section [iVCl 



C. 1 + 2 — > 4, endpoint contribution 

Finally, there is a possibility that both partons emerg- 
ing from the perturbative splitting in mD experience 
hard collisions straight away, without any further evo- 
lution. Jet imbalances stemming from such an eventual- 
ity are no longer independent but, on the contrary, are 
strongly correlated. At the "Born level" one has 
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With account of additional radiation, the delta-function 
gets replaced by the pole enhancement of the differential 
cross section in the imbalance of imbalances in a spe- 
cific region of jet momenta when the pair imbalances are 
practically equal-and-opposite: : 



da 



<P5, a d?6- 



oc 



2 1 



d 2 S' 2 ' 



«5 2 



$13 



Importantly, this contribution is also double-collinear en- 
hanced and therefore has to be taken into full consider- 
ation. This eventuality is not incorporated in Eq. (26) 
and should be taken care of separately. 

In this kinematical region the parton that carries the 
compensating transverse momentum —5' can be pro- 
duced off one of the external legs, as shown in Fig. |4| 




The corresponding contributions to the differential 4-jet production cross section can be combined into the following 
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relatively compact expression: 



da 



n 2 da { 2 3 ^ 

d 2 Sis d 2 S 2 4: dii dii 



part 



P. 



1,2 



a s (6 2 

2tt S 2 ' \ x± + X 2 
/2 .G c a (x 1+ x 2 ;S' 2 ,Q 2 ) 



Si(Q 2 ,S 2 )S 2 (Q 2 ,S 2 ) 



Xi + X-2 



S 3 « 2 ,i' 2 )54W 2 ,i' 2 )x[ 2t Dffe^4;{' 2 ,f 2 )| + {a^b}. 
I 



(27) 



Once again, integral over A on the r.h.s. is implicit here. 
In this formula the splitting of the parton "c" from the 
hadron "a" is written explicitly, represented by the split- 
ting function P c 1,2 . The splitting occurs at the virtu- 
ally scale S 2 . However, accompanying production of 
secondary real particles is vetoed starting from a much 
smaller scale, 8' 2 -C S 2 . This explains appearance of the 
Sudakov form factors of five participating partons in ( 27 1 . 



In the impulse approximation, there is one parton pro- 
duced with the transverse momentum —6' compensat- 
ing the 4-jet imbalance, while transverse momenta of all 
other real partons are smaller: <C S' 2 . Production 
probability of this parton integrated over its energy (ra- 
pidity) is embedded into the derivative of the product of 
three Sudakov form factors and of the initial parton dis- 
tributions depending on 5' 2 as the upper evolution scale. 



As before in Eq. (25), by differentiating form factors we 
obtain soft gluon radiation off three external lines "c" 
("0"), "3" and "4". Differentiation of the parton distri- 
bution functions describes real parton production due to 
"hard" splittings. 



The full answer is given by the sum of Eqs. (25)-(27l 



7T 2 da (i 



4) 

d 2 Si 3 d 2 8 2i 



(3->4) 

l 

d 2 5iz d 2 5 2i 



tt 2 da\ 



(3-C4) 



7T 2 da.^ 
d 2 d 13 d 2 S 2 



(28) 



D. Distribution over jet imbalances in the extreme 
back-to-back limit 8 2 ,8' 2 -> 



As we have repeatedly stressed above, double parton 
collisions constitute but a small correction to 4-jet pro- 
duction if characterized by their contribution to the total 
cross section. However, "higher twist" four-to-four and 
three-to-four processes are specific in that they populate 
the kinematical region of relatively small pair jet imbal- 
ances, S 2 3 ,S 2 4 , S' 2 = (<5i3 + S 24 ) 2 < Q 2 - 

Being a multi-scale hard process, back-to-back jet 
production cross section possesses double logarithmic 
parton form factors. The Sudakov suppression is the 
price one pays for having the differential distribution 
peaked, da oc S~ 2 , which enhancement implies vetoing 
the bremsstrahlung radiation of gluons with transverse 
momenta in a broad interval ranging from k\ > S 2 all 
the way up to k\ < Q 2 . 



The formulas ( 25 )— ( 27 ) for the differential imbalance 



spectra were derived in the impulse approximation in 
which it is a single parton (soft gluon) that compen- 



sates the transverse momentum imbalance, while radi- 
ation of other partons with k\ > 5 2 is strictly vetoed. 
However, in the extreme S 2 — > limit the form factor 
suppression becomes so strong as to override the pole 
enhancement. In these conditions the impulse approxi- 
mation is no longer valid, and instead of vetoing accom- 
panying radiation one has to impose the condition that 
the real partons carry, in aggregate, the given transverse 
momentum: fc_u = —5. This implies exponentiating 
soft gluon radiation in the impact parameter space, and 
evaluating the inverse Fourier transform to access the re- 
sulting 5-spectrum [33ti35j . The result is the flattening 
of the spectrum in the S 2 — > limit. 

The precise width of this plateau depends on par- 
ton form factors as well as on parton distributions, and 
should be calculated numerically. However, a rough semi- 
quantitative estimate of the value of the critical trans- 
verse momentum p where the spectrum starts to flatten 
out can be obtained from the following simplified consid- 
eration. 

In the double logarithmic approximation one can ne- 
glect effects due to parton distributions and state that the 
plateau develops when the product of relevant Sudakov 
form factors just compensates the pole enhancement: 



din / 



From Eqs. @ we obtain the condition 

^ « s (Po) 



In— = 1, 



7T p 



(29) 



where £c is the sum of the "squared color charges" 
(Casimir operators) of participating partons. Substitut- 
ing the one-loop expression for the running coupling, we 
have 



hi 



In 



Q 2 



A 2 A2 
Jl QCD ^QCD 



2S r 



(30) 



with f3 2 = HN c /3 — 2nf/3 the leading coefficient of the 
QCD /3-function. The estimate of the "flattening mo- 
mentum" follows: 



2 A 2 ( QM (, , fh 



(31) 



Using /3 2 = 9 for nj = 3 light quark flavors, the exponent 
7 in Eq. (31) equals 0.372 for a quark pair (Sc = 8/3), 
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0.456 for a system of a quark and a gluon (So = 13/3), 
and 0.571 for two gluons (Xc = 6). 

The value po characterizes a typical total transverse 
momentum of the final state produced in a hard 2-parton 
collision, be it a Drell-Yan pair, W/Z or a system of four 
hadron jets (or, say, three jets and a large transverse 
momentum photon). 

In the case of four- (three-) parton collisions, one has 
double back-to-back enhancement, accompanied by the 
four Sudakov form factors, in place of two. So, the pat- 
tern of the spectrum flattening in the S 2 ,5 12 — > limit 
remains roughly the same. 

We conclude the discussion of the differential distribu- 
tion with two remarks. 

As we have seen above, in the back-to-back kinematics 
the 4-jet production in double parton collisions (4 — » 4, 
3 — > 4), and due to higher order QCD effects in stan- 
dard two-parton collisions (2 — > 4), are comparable. This 
makes the experimental studies of double hard collisions 
not an easy task. To extract MPI, one should learn to 
reliably subtract the 2 — > 4 contributions. This "back- 
ground" should be theoretically predicted and generated 
at the NLO level (for the amplitude) in order to accom- 
modate loop correction effects discussed above in Sec- 
tion hvb 

On the experimental side, the Tevatron experiments 
[TBHTS] have set a fashion of searching for MPI by look- 



ing at angular correlations between transverse momen- 
tum imbalances in 7 + 3 jets (and 7 + 2 jets, [18]) events. 
Such a strategy can well be used to signal the presence 
of MPI. At the same time, it is not suited for quantita- 
tive analysis as it mixes together contributions of hard 
and soft physics. In order to extract and study dou- 
ble hard collisions, and thus to get hold of inter-parton 
correlations inside hadrons, one should instead base the 
measurements on the values (and relative geometry) of 
transverse momentum imbalances. 



VI. TOTAL CROSS SECTION 



Integrating the differential distribution ( 28 ) over im- 
balances d 2 8i3, d 2 S24 we obtain the total cross section of 
the two dijet production in the double parton collision 
process in the leading collinear approximation: 



da{x\, X2, £3, X4) da 3 da 



T 24 



dti dt2 



dii dt2 



1 

5*4 



1 



(32a) 



Here £4 and S3 are Xi- and Q 2 -dependent parameters 
that describe effective interaction areas characterizing 4- 
and 3-parton collisions, correspondingly. They are given 
by the following expressions: 



S 4 1 (xi,x 2 ,x 3 ,X4;Q 2 ) 



d 2 A 



2]D a {x\ 1 x2\Q ,Q , A) [2 ]D 6 (x 3 , x 4 ; Q , Q ; -A) 



[2}D a {xi,x 2 ]Q ,Q ;A) [1 ]D b (x 3 ,x 4: ;Q ,Q ; -A) + [ 1 ]D a (xi,x 2 ;Q , Q ; A) [ 2 pb(%3, £4; Q ,Q ;-A) 



(32b) 



and 



S 3 1 (x 1 ,x 2 , 



X 3 ,X4 



^ G c a (x 1 +x 2 ,6 2 ,Q 2 ) n3 ,4 



(32c) 



X\ + x 2 



[2 p^(x 3 ,x A ;S 2 ,S^A) + (a 06; l,2o 3, 4) 



Unlike the case of 4 — > 4, at the level of the differen- 
tial spectrum the contribution 3 — > 4 due to an interplay 
of non-perturbative and perturbative two-parton corre- 
lations is not given by an expression with full derivatives 
over two relevant scales; therefore the final expression for 
SV contains logarithmic integration over the hard scale 
5 2 running up to the overall hardness scale of the parton 
scattering Q 2 . 



VII. CONCLUSIONS AND OUTLOOK 

In this work we aimed at generalization of the QCD 
factorization theorem to multiple parton interactions in 
hadron-hadron collisions. Such generalization is neces- 
sary for carrying out systematic studies — both theo- 
retical and experimental — of the parton correlations 
inside hadrons, and allows one to incorporate in a model 
independent way higher order pQCD effects. In the 
small- cc domain, generalization of the MPI analysis pre- 
sented here may be looked upon as realization of the 
Gribov Pomeron Calculus, with exact account of energy- 
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momentum conservation. 

We considered four jet production and have shown that 
double hard parton collisions play a significant role in the 
back-to-back kinematics. This kinematics one selects by 
demanding the transverse momentum imbalances of two 
pairs of jets, 5 2 3 = (j ±1 + j± 3 ) 2 and <5f 4 = {j± 2 +J±4) 2 , 
to be much smaller than the overall hard scale of the 
scattering process: ^13,^24 

We were able to analyze and calculate corresponding 
differential spectra and total cross sections with account 
of all collinear enhanced pQCD radiative effects giving 
rise to Sudakov form factors and to scaling violations in 
one-parton and generalized double parton distributions, 
2GPDS, that have been introduced in pQ. Our main re- 
sults are Eqs. (25)-(27) for the differential distribution 



and Eq. (32) for the total cross section of 4-jet produc- 



tion in the back-to-back kinematics. They are derived in 
the leading logarithmic (collinear) approximation. The 
formalism developed here provides the framework, and 
a starting point, for systematic calculation of NLL and 
NLO corrections that may be needed for more careful 
account of the MPI effects at LHC. 

In particular, in the present study we did not hunt 
for single logarithmic corrections due to large-angle soft 
gluon radiation. These (formally subleading but poten- 
tially important) logarithmic effects depend on details of 
the color transfers in hard parton scattering processes. 
They should be treated together with a delicate question 
of algorithmic determination of transverse momenta of 
the hadron jets. The corresponding study [UJ exists for 
the single hard scattering case, 2 — > 2, and provides a 
solid base for analogous treatment of large-angle gluon 
radiation in double parton collision processes. 

We demonstrated that, in addition to the QCD- 
improved parton model picture of 4-parton collisions, 
pQCD also reserves an important place to 3-parton colli- 
sions producing four jets. Such 3 — > 4 transitions are the 
novel hard double scattering processes that are system- 
atically treated in this work for the first time. 

3 — > 4 parton process appears when a parton from 
one colliding hadron splits into two partons that enter 
into double hard interaction with two partons from the 
(non-perturbative) wave function of the second hadron. 
Perturbative parton splitting occurs at relatively small 
impact parameters, while the relative impact parameter 
distance between two hard vertices in a 4 — > 4 collision 
can be as large as the transverse size of the hadron. Para- 
metrically, 3 — > 4 and 4 — > 4 subprocesses are compara- 
ble (they are of the same order in a s and have similar 
structure of the back-to-back enhancement). At the same 
time, the difference in geometry of the interaction pro- 
cesses results in a numerical enhancement of 3 — > 4 with 
respect to 4 — > 4 subprocess (by a numerical factor of 
4-5). As a result, 3 — > 4 transitions should contribute 
significantly to 4-jet production, even if they happen to 
constitute formally a small correction to a 2GPD. 

It is important to stress that the processes where 
one parton splits perturbatively into two in both col- 



liding hadrons do not contribute to the 2-dijet pro- 
duction, since the corresponding differential distribution 
lacks the necessary back-to-back double pole enhance- 
ment da/(d 2 5) 2 oc &~ A . In fact, such configurations rep- 
resent a one- loop correction to the 2 — > 4 jet produc- 
tion mechanism and as such belong to the class of the 
standard single hard interaction processes rather than 
to MPI. This result of [T] has been later confirmed in 
Refs. [H1IMI- 

The role of 3 — > 4 parton subprocesses is twofold. 

When the perturbative splitting occurs at momentum 
scales smaller than the pair jet imbalances 8 2 (that is, 
"deep inside the DGLAP ladder"), such contributions 
exhibit the same two-pole enhanced structure as the con- 
ventional 4 — > 4 transitions already discussed in pQ : 



da 



1 1 



d 2 <5i3 d 2 S 2 4 613 6% 4 



The corresponding differential distribution is a direct 
generalization of the "DDT formula" 20J containing 
derivatives over transverse momentum imbalances, with 
Eq. p5| describing 4 — > 4 collisions and Eq. (26) — a 



crosstalk between large- and small-distance two-parton 
correlators. 

At the same time, the parton splitting may occur just 
at the scale S 2 and produce two partons that, without 
any further evolution, interact with two partons from the 
wave function of the second hadron. Such eventuality 
also has a necessary back-to-back enhancement but of a 
different structure: 



da 



1 1 



d 2 S 13 d 2 S 2 4 S' 2 



=1 



24- 



Such contributions are absent in the parton model. They 
cannot be represented as the product of two 2GPDS and 



must be taken into account separately, see Eq. ( 27 ). A de- 
tailed numerical investigation of these contributions will 
be presented elsewhere [3"B"] . 

Formulas derived in this work can be used to address 
various aspects of the LHC physics. Our results are di- 
rectly applicable to the case of double Drell-Yan pro- 
cesses at LHC. In particular they can be applied to the 
recently observed production of heavy quarkonia at the 
LHC [37J . The formalism developed in the present paper 
should help to improve the accuracy of the prediction of 
QCD backgrounds in searches for new physics. It also 
can be used for further studies of the nucleon structure, 
in particular of the short-range non-perturbative inter- 
parton correlations in the nucleon and the nuclei, for the 
study of the higher twist contributions that are unavail- 
able in the conventional DIS processes, etc. In addition, 
it would be interesting to extend the formalism devel- 
oped here to investigate QCD medium effects that mani- 
fest themselves in AA collisions in recent experiments at 
CERN. 

On the experimental side, observation of the parton 
splitting processes discussed in this paper will give the 
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first direct evidence of an interplay between large- and 
short-distance QCD correlations in hard processes. To 
achieve this goal one has to look for enhancement of 4- 
jet production in the kinematical region of jet transverse 
momenta 

(ji± +J2X + Ja_L + 34±) 2 < (ji± + J3_l) 2 - (.h±+j4±) 2 - 
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2. Amplitudes 



Adding the virtual propagators 

a x 2 , x\ 

W sf 4 and Ws' 

correspondingly, and extracting the factor (xi + x 2 + 
q)(xi + X2) we get the amplitudes 

n t aS' x 2 S 13 - XxS 24 , . 

= T • — x P (A.3a) 

o u xi + x 2 

1= Tl . ^ ~ ° g 13 x £2^24 . p (A . 3b) 



X\ + x 2 + a 



n _ rr, X 2^' + aS 24 Xl(5l3 

2 — 1 2 ' , , X ~+2 ' r - 

xi + x 2 + a of 3 



(A.3c) 



Appendix 

Here we present the analysis of real parton production 
in the course of 1 — > 2 parton splitting, — >• 1 + 2, de- 
termining the correlated 4-jet production. We consider 
three contributing diagrams: that of Fig. ga) for emis- 
sion off the initial parton "0", and those of Fig. [3] for 
emission off the offspring partons "1" and "2". 



1. Vertices 

It is straightforward to write down exact expressions 
for the amplitude using the economic technique of quasi- 
parton states developed by Bukhvostov, Frolov, Lipatov 
and Kuraev (BFLK; not to be confused with BFKL) and 
described in [35]. In the BFLK technique the quark and 
(axial gauge) gluon Green functions are replaced by prop- 
agators of two physical ± helicity states. In general, in 
this framework additional terms looking like contact four- 
parton interaction vertices arise. They, however, do not 
produce collinear enhanced contributions. 

Effective three-parton interaction vertices are propor- 
tional to the linear combination of the transverse mo- 
menta of participating partons weighted with the longi- 
tudinal light-cone momentum fractions: 



(A.l) 



In the matrix element the vector index is contracted with 
the gluon polarization vector (for details see [38]). 
The products of two vertices for our amplitudes read 

: To ■ (xi + x 2 + a)5' x (x 2 6 13 - x^m) (A.2a) 

1 : Ti • (-(xi+a)<5i 3 - xrfu) x (xi+x 2 )5 2 4 (A.2b) 

2 : T 2 ■ ((x 2 + a)8 2i + x 2 <5 13 ) x (x 1 +x 2 )5 13 (A.2c) 

Here the first vector corresponds to the real parton (soft 
gluon) radiation, the second vertex — to the internal 
hard splitting. 



Here P is the common propagator factor defined in ( 23 1 



P- 1 =x 1 <5 2 2 4 +x 2 ^3 + (x 1 +x 2 )r- 1 ( 5' 2 , r^ °^ c l± x ^ . 

X\X 2 

(A.4) 



3. Limits 

In the case of soft gluon radiation, a -C Xi, the pa- 
rameter r defined in (A.4) becomes numerically small. 



Therefore one should examine, separately, two kinemat- 
ical regions: 

1. 6' 2 « r ■ S 2 

Only "0" contributes: 



M ~ To 



2. S 2 > 5 12 ^>r-6 2 



>13 



aS' 



S' 2 (X! + x 2 )S 2 3 



(A.5a) 



Here the amplitudes (A.3b) and (A.3c) contribute, and 
we obtain 



M^(T 1+ T 2 ).^^^- 
x\ + x 2 6{ 3 x\x 2 o' 1 



(A.5b) 



reproducing the structure of (A.5a). Emissions from 1 



and 2 are coherent, and due to conservation of the color 
current, T\ + T 2 = To, the two expressions coincide. As 



a result, (A.5a) applies in the entire kinematical region 
S' 2 <C S 2 , irrespectively to the value of r. 

This contribution — real parton production off the in- 



coming line"0" — isapart ofthc3 — > 4process, Eq. (27l, 
described in Section IVBI 



Additional collinear enhanced contributions arise in 
two complementary regions. 

3. 5 2 4 < <5 2 3 ~ S' 2 

Here the amplitude 1 dominates: 



M 



S 13 S\ 



13 024 



x\ + x 2 + a S 2 3 S24 



(A.6a) 
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4. 6 2 13 « 51 



Alternatively, here the dominant amplitude is that with 
emission off the line 2: 



M 



T 7 - 



a 



xi + X2 + a Sf 3 5? 



$13 $24 
13 °24 



(A.6b) 



The contributions (A.6| are related with internal evolu- 



tion of the 2GPD of the hadron a and are contained in 
Eq. (25). Specifically, in the kinematical region where 



the imbalances of two jet pairs are significantly differ- 
ent, these contributions correspond to the case when the 
differentiation over the smallest of the two imbalances 
applies to the 2GPD and acts upon the integration limit 
of the correlation term in (18). 
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